Fractional decay of quantum dots in real photonic crystals 
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We show that fractional decay may be observable in experiments using quantum dots and photonic crystals 
with parameters that are currently achievable. We focus on the case of inverse opal photonic crystals and 
locate the position in the crystal where the effect is most pronounced. Furthermore, we quantify the influence 
of absorptive loss and show that it is a limiting but not prohibitive effect. 
OCIS codes: 270.1670, 270.5580. 



Spontaneous emission is a resonant process in the sense 
that an emitter interacts with modes of the electromag- 
netic field spectrally close to the electronic transition 
frequency. Moving the emitter to another medium or lo- 
cation at which the field strength of the electromagnetic 
vacuum modes differ will lead to changes in the sponta- 
neous emission. In most cases the light-matter coupling 
is weak and the emitter decays exponentially in time 
with different decay rates at different locations. This is 
the weak coupling Purccll regime. 

In certain cases the coherent coupling of the emitter 
to a highly structured electromagnetic vacuum leads to 
non exponential decays. In particular, a regime of so- 
called fractional decay has been pointed out [1,2]. In 
this regime, the emitter coherently interacts with modes 
of low group velocity in such a way that it never fully 
decays, but rather remains in a superposition of the ex- 
cited state and the ground state. This may happen in 
media with rapid variations in the spectral and spatial 
distribution of electromagnetic modes, as described by 
the local optical density of states (LDOS) [3]. Photonic 
crystals offer the ability to manipulate the LDOS and 
change it as compared to the case of a homogeneous 
medium. The Purcell effect has been shown experimen- 
tally for a variety of emitters, e.g, quantum dots (QDs) 
in inverse opal photonic crystals [4] , but there is to date 
no demonstration of fractional decay. 

In this Letter we intoduce a practical measure of the 
degree of fractional decay and use it to investigate de- 
cay dynamics of QDs near the band edge of a photonic 
crystal. For calculations of decay dynamics we follow the 
approach of Vats et. al [5] . Contrary to the general treat- 
ment in [5] we focus in this Letter on the possible re- 
alization of fractional decay using specific and realistic 
structures. In particular the investigations are based on 
the actual LDOS of a three dimensional photonic crystal 
obtained from plane wave calculations and extended to 
include also effects of absorptive losses, see Fig. 1. Ab- 
sorption is shown to be a limiting factor and we present 
quantitative results showing the degreee of fractional de- 
cay achievable for available QDs and practically relevant 
material loss. 
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Fig. 1. Projected LDOS at the H point of a closed packed 
Si inverse opal in units of Wg ff /(37r 2 c 3 ) calculated with 
plane wave expansion (PW) and analytical approxima- 
tions in the vicinity of the band edge for two different 
losses. Inset shows the unit cell of an inverse opal and 
high symmetry points. 



The QD is modelled as an initially excited two-level 
system with transition frequency oj eg and we write the 
general state of the coupled electron-photon system as 
|*) = c e (t) |e,0) +X^ c <?,m (*)lff>M)> where |e,0) denotes 
the electron in the excited state and no photons and 
\g,jj) the electron in the ground state and one pho- 
ton in mode /z. The time evolution is governed by the 
Schrodingcr equation, iH\^f) = H \^>), with the Hamilto- 
nian H in the dipole and rotating wave approximation. 
Using a Laplace transform [5] the equations of motion 
are solved in the frequency domain to yield 



c e (r,w) 
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where r is the QD position and ui = to/ oj eg is the scaled 
frequency. We have split the spectrum into a pole term 
with residual a_i, pole position ujq and a rest term. The 
dimensionless parameter (3 = T$/ (2iruj eg ) is the vacuum 
decay rate, To, scaled by the transition frequency, and is 
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given as (3 = q 2 p 2 / (6hm 2 ir 2 eo c 3 ), where q,p,h,m,eo 
and c denote electron charge, momentum matrix ele- 
ment, reduced Planck constant, electron mass, free space 
permittivity and vacuum speed of light, respectively. Ex- 
perimental values range from (3 w 10~ 8 for InAs QDs [6] 
to (3 w 6 x 10~ 8 for PbSe QDs [7] with so-called inter- 
face defect QDs possibly reaching values of (3 ~ 10 -6 [8]. 
The function G(r, uj) is given for frequencies above the 
integration path in the complex plane as 



G(r,(D) = iuj 



x 2 (uj — x) 



(2) 



in which p p (r,uj) = p p (r,uj)/po(uj e g) is the ratio of the 
LDOS to the vacuum LDOS at the emitter frequency. 
The decay is a resonant process and is governed by the 
LDOS in a narrow frequency interval around uj = 1. 
The remaining LDOS, however, does contribute an over- 
all Lamb shift of the spectrum. To model this effect we 
include an integration using p p = po for uj < 0.95 and 
1.01 < uj. A cutoff is chosen at ujc = 10 5 , correspond- 
ing to the Compton frequency [5]. For 0.95 < w < 1.01 
the integral is carried out using the accurate LDOS 
Pp = Pbe, as calculated below. 
The projected LDOS is defined as 



P P (r,w) 



E. 



lep-E^r)! 2 ^-^), 
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where the sum is over all modes of the electromagnetic 
field indexed by p and e p is the orientation of the emit- 
ter. The functions E M (r) = (r|E M ) denote the spatial and 
spectral distribution of the modes and are normalized as 
(E a |efl( r )|E @)v = $a,i3, where €r is the relative permit- 
tivity and V is the normalization volume. In vacuum the 
LDOS, Eq. (3), is given as Po {uj) = uj 2 /(3tt 2 c 3 ). 

Fig. 1 shows a zoom in on the LDOS, p p (uj), at the 
H point of an Si inverse opal (en = 11.76) close to the 
upper edge of the band gap. The LDOS was calculated 
using a method similar to that of [9], corrected for the 
reduced symmetry of the electric field [10] and using 169 
plane waves and 1232944 k-points distributed over half 
the full Brillouin zone. The plane wave approach results 
in a discrete sampling of the LDOS, which effectively 
limits the slope of the sampled LDOS and leads to in- 
correct results when used in calculations of fractional 
decay. For this reason, and in order to include losses in a 
pcrturbative way, we analyze analytically the LDOS in 
the vicinity of the band edge. 

The upper band edge in inverse opals is defined by the 
ninth band at the so-called X point only [9] . Considering 
the contribution from just a single band, we follow [11] 
and rewrite the sum, Eq. (3), as 



P P (r,u) 



V 



(2tt)= 



E k ,,(r) 



|Vw(k)| 



■dk, (4) 



in which the sum is over two different polarizations and 
the integration is over the dispersion surface of constant 
frequency uj corresponding to the ninth band only. We 



now expand the integrand in Eq. (4) in powers of k and 
carry out the integration. To lowest order the LDOS is 
given as p BE (u) = K BE (r)y/uj - uj be , where uj be is 
the band edge frequency and the band edge parame- 
ter, K BE (r), is related to the curvature of the disper- 
sion surfaces and the projected electric field at the point 
r. The black curve in Fig. 1 illustrates how the square 
root is indeed the limiting form of the LDOS close to 
the band edge. Fig. 2 shows values of K BE along lines 
between symmetry points of the Wigner-Seitz cell of a 
nearly closed packed Si inverse opal (hole radius per lat- 
tice constant, R/a = 0.3436). The analytical approach 
allows for the use of only 5 k— points in each direction 
for the determination of the curvature and 1243 plane 
waves to achieve convergence [12]. 
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Fig. 2. The parameter K BE (in units of po(ui eg )/uj e g ) 
as a function of position in the Wigner-Seitz cell (shown 
in the inset) for the three principal emitter orientations. 
Grey dash-dotted line shows values of K BE at the H 
point for different R/a (top scale). 

Introducing loss in the material, e B — > fjj + if/. leads 
to a shift in frequency. For small losses we use first order 
perturbation theory [13] to write uj = uj 1 -- ^ — iS, where 
is the frequency in the absence of losses and 
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where subsript C denotes the volume of the lossy mate- 
rial only leading to / = (E M |e fl |E A( )c/(E M |e iJ |E M )v < 1. 
For nonzero 5 we rewrite the band edge LDOS as [11] 



p BE (r,uj) = K BE (r) / 

J UJl 



y/x - UJ BE 



5/tt 



(uj - x) 2 + S 2 



■ dx, 



which shows that the effect of absorption is to broaden 
the modes as well as to introduce states below the upper 
edge of the band gap (dashed curve in Fig. 1). 

Using the above expression for p BE (r,uj), the spec- 
trum, c e (r, uj) is calculated from Eqs. (1) and (2). The 
temporal evolution is subsequently obtained by trans- 
formation back to the time domain. Under this transfor- 
mation, the pole term in Eq. (1) is conveniently handled 
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analytically, resulting in a decreasing exponential part. 
The absolute square of the residual denotes the strength 
of the pole term and is equal to the value of the expo- 
nential part at t = 0. The case of |a_i| 2 = 1 results in 
c'(lo) = and the spectrum consists of only a single pole 
term. This is characteristic of the weak coupling Purcell 
regime and the decay is exponential with a decay rate 
r = Top p (uj e g) / po(ui e g). On the other hand, |a_i| 2 < 1 
results in a non-zero rest term and consequently a devi- 
ation from the Purcell regime and we define |o_i | 2 < 1 as 
the condition for fractional decay. The residual depends 
critically on the light-matter coupling strength, (3Kbe, 
relative to the absorption. The former depends on the 
QD as well as position in the photonic crystal, cf. Fig. 2. 

For a specific example we consider now colloidal PbSe 
QDs, emitting at v PbSe w 1.3 ■ lO^s" 1 (fi w 5.5 x 1(T 8 
[7]) and placed at the H point in a closed packed Si 
inverse opal [Kbe ~ 10). Fig. 3 shows the resulting de- 
cay curves for different absorption. For vanishing losses, 
the population tends to a non-zero value at long times 
with |a_i| 2 = 0.84. At small finite losses a fractional ef- 
fect is still visible with |a_i| 2 = 0.87 at an absorption 
length of a = 3 x lO^cm -1 (S/uj PbS e = 10" 10 ) and 
|a_i| 2 = 0.96 at a = 3 x lO^cm -1 {5/u PbSe = 1CT 9 ). 
We note that absorption in Si at this frequency may be 
as low as a w lO^cm" 1 {5/uj PbSe w 10" 13 ) [14]. Fig. 
3 therefore shows that fractional decay is observable for 
real QDs in dielectric photonic crystals exhibiting ab- 
sorptive losses. 
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Fig. 3. Decay from PbSe QDs at the H point in an 
Si inverse opal (/3K P e = 5.5 x 10 -7 ) at detuning 
wss/wpbSe = 1 — 8.309 x 10~ 6 and different losses. 
Dashed curves show exponential parts only. 

For a given system the detuning of the emitter rela- 
tive to the band edge defines the exact modes, and hence 
group velocity, of the emitted light. Therefore, the resid- 
ual depends also on the detuning and we define the pa- 
rameter Df as the minimum value of |a_i | 2 for optimized 
detuning. In Fig. 4 we show Df as a function of (3K P e- 
The curves were obtained for each f3K p e by varying the 
detuning until a minimum was found. The figure shows 
that a profound degree of fractional decay is possible for 
a range of experimentally relevant material parameters. 

In conclusion we have used an analytic expression to 
the band edge LDOS to investigate fractional decay dy- 
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Fig. 4. Degree of fractional decay obtainable versus light- 
matter coupling strength, (3Kbe, at different absorption. 

namics in inverse opals. The analysis has revealed the 
position in the crystal that is most suitable for observa- 
tion of fractional decay. Furthermore, we have extended 
the analysis to include absorptive losses and calculated 
the degree of fractional decay obtainable for given losses 
and light-matter coupling strengths. The analysis shows 
that absoption has a limiting but not prohibitive effect 
and that fractional decay may be possible to achieve us- 
ing, e.g., PbSe QDs in Si inverse opals. 

Enlightening discussions with Willem Vos and Niels 
Asgcr Mortensen are greatly appreciated. 
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